SCIENTIFIC 

REPpRTS 




OPEN 



SUBJECT AREAS: 

ELECTRONIC PROPERTIES 
AND MATERIALS 

NONLINEAR PHENOMENA 



Received 
19 November 2013 

Accepted 
17 March 2014 

Published 
2 April 2014 



Correspondence and 
requests for materials 
should be addressed to 
J.Z. (zhoujunzhou@ 
tongji.edu.cn) 



Wave-packet rectification in nonlinear 
electronic systems: A tunable 
Aharonov-Bohm diode 



Yunyun Li\ Jun Zhou\ Fabio Marchesoni^-^ & Baowen Li^ - 



^Center for Phononics and Thermal Energy Science, School of Physics Science and Engineering, Tongji University, Shanghai 
200092, People's Republic of China, "^Department of Physics, University of Camerino, 1-62032 Camerino, Italy, ^Department of 
Physics and Center for Computational Science and Engineering, National University of Singapore, Singapore 1 1 7456, Republic of 
Singapore. 

Rectification of electron wave-packets propagating along a quasi-one dimensional chain is commonly 
achieved via the simultaneous action of nonlinearity and longitudinal asymmetry, both confined to a limited 
portion of the chain termed wave diode. However, it is conceivable that, in the presence of an external 
magnetic field, spatial asymmetry perpendicular to the direction of propagation suffices to ensure 
rectification. This is the case of a nonlinear ring-shaped lattice with different upper and lower halves (diode), 
which is attached to two elastic chains (leads). The resulting device is mirror symmetric with respect to the 
ring vertical axis, but mirror asymmetric with respect to the chain direction. Wave propagation along the 
two diode paths can be modeled for simplicity by a discrete Schrodinger equation with cubic nonlinearities. 
Numerical simulations demonstrate that, thanks to the Aharonov-Bohm effect, such a diode can be operated 
by tuning the magnetic flux across the ring. 



The design of current rectifiers, also known as diodes, has spawned widespread interest in many branches of 
physics. The analogs of electronic diodes, such as thermal diodes^ ^ optical diode^'^ and acoustic diodes^'^, 
have been theoretically proposed to effectively manipulate energy transport and wave propagation. Thanks 
to recent advances in nanotechnology, prototypes of some of them have indeed been realized^°'^\ The operation of 
this class of quantum devices requires at least two ingredients, namely, nonlinearity and spatial asymmetry along 
the direction of propagation^^ (longitudinal asymmetry). A similar conclusion was reached in the classical theory 
of Feynman's ratchets and Brownian motors, where autonomous transport of neutral particles and other small 
objects in nonlinear asymmetric devices is sustained by ambient fluctuations^^. 

In electronic systems the interaction between electrons and lattice vibrations, or phonons, is responsible for 
appreciable nonlinear transport effects^^"^^. This scenario applies, for instance, to conductive polymers, where 
solitons play a fundamental role in the charge-transfer doping mechanism^^; to macromolecules such as DNA, 
where the nonlinear electron-phonon coupling is likely to play an important role^^"^\ In this context the dynamics 
of an electron wave-packet subject to nonlinear electron-phonon interactions was conveniently modeled in terms 
of discrete nonlinear Schrodinger equations (DNLSE)^^. The DNLSE's have also been employed to describe the 
propagation of longitudinal Block waves in layered or phononic crystals^^ light through networks of coupled 
optical waveguides in Kerr media^^ and interacting cold atoms on an optical lattice^^. 

A composite quasi- one dimensional nonlinear device inserted along a transmission line operates as wave diode 
when its transmission coefficients depend on the direction of the incoming wave. The required mechanism of 
longitudinal symmetry breaking is typically established by grading the composition and geometry of the diode 
components^^'^^. Moreover, for practical purposes any such device must be tunable, which raises the question of 
how to control the diode output by means of an external signal. 

In this paper we tackle both issues of asymmetry and control, by considering a novel diode topology. A sketch of 
the device we propose is shown in Figure 1, where a ring (diode) is connected to two chains (leads) parallel to one 
of its diameters. Both the chains and the ring are modeled respectively as elastic and nonlinear discrete lattices. 
The upper and lower halves of the ring have different propagation properties. With this design the left-right 
mirror (or longitudinal) symmetry of the device is preserved, whereas its upside-down (or transverse) symmetry 
is broken. Under the operating conditions discussed in the earlier literature, such a device would be of no use as a 
wave rectifier^^. However, this is no longer true after a magnetic field perpendicular to the ring plane is switched 
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Figure 1 | Sketch of an AB wave diode: an upside-down asymmetric, nonhnear ring-shaped lattice is connected to two elastic chains with Np = 500 and 

Ni = 2. The magnetic field is applied perpendicularly to the ring. The ring site interactions have parameters s+ = So(l ± Se) and g+ = go(l ± 6g), 
see equation (1), where ± denote the upper (lower) half ring, 5g and 5g control their asymmetry, and go sets the degree of nonlinearity in the ring. At the 
junction sites, n = Ni, 6e = 5g = 0. 



on. Owing to the nonlinear nature of the ring, the difference between 
the Aharonov-Bohm (AB) phases electrons acquire along the upper 
and lower paths depends on the side they enter the diode. We prove 
that, as a consequence, the left- right and right-left transmission coef- 
ficients across the AB ring can be different and the relevant rectifica- 
tion factor can be easily controlled by tuning the magnetic field. 

Results 

Model. The nonlinear propagation properties of the ring were 
modeled by the following DNSLE for the electron wave function \\f^^^ 

i^n=^n^n+gn\^nHn-^in + m,n^n + m- (l) 
m 

The sites of the left (right) chain are labeled from — Np to — Ni — 1 
(from Ni + 1 to Np); the sites of the upper (lower) ring half are labeled 
from (-Ni + 1) to (Ni - 1) [from (-Ni + 1)' to (Ni - 1)'], and the 
ring-chain junction sites are ±Ni. Sn represents the on-site energy 
and the constant gn is meant to reproduce the local nonlinear 
electron-phonon coupling^^. In particular we set: Sn = gn = 0 on 
the elastic chains; Sn — £± — £o(l — Se)^ gn — g± — go(l — ^g) in the 
upper (lower) ring sides; and, finally, £+]V/ = eo and g±Ni=gQ at 
the junction sites. The parameters 5g and 5g were used to control 
the asymmetry between the upper and lower halves of the ring, while 
go was tuned to vary its overall nonlinearity. 

In the last term of equation (1) the summation over m was 
restricted to hopping between nearest neighbors, for which tn+m,n 
= 1. All sites are coupled to two nearest neighbors except the junc- 
tion sites with n = ±Ni, which have three nearest neighbors with 
n + m = ±(Ni - 1), ±(Ni 1)', and ±(Ni + 1). The hopping rates 
between all other sites were set to zero, so that there is no direct 
interaction between the upper and the lower side of the ring. 

The effect of the magnetic field on a charge q hopping along the AB 
ring was accounted for by means of the so-called "Peierls substi- 
tution"'°''\ 

W.,n-^n + .,nexp(^±/^^J, (2) 

where Nr = 4Ni is the number of ring sites, O the magnetic flux 
connected to the ring, = h/q the quantum of magnetic flux and ± 
stand for the positive or negative orientation of the n ^ n + m hop. 
The phase modulation in equation (2) is restricted to the hopping 
rates between ring sites. 

We further notice that for our choice of tn+m,n equation (1) 
can be easily derived from the Hamiltonian = 
^{y[^n]-i^n + \Vn + Vn+i^n))^ where the quartic on-site poten- 

n 

tial = &n\^nY I 'An 1^/2 is nonzero only inside the ring. 

Accordingly, both the energy and the norm of the propagating 
wave-packet are conserved, as expected for electrons propagating 
in a non-dissipative medium. In our simulations we focused on the 
case of incident Gaussian wave-packets. 



^„(0) = /oexp(^-^^^^^-*onj, (3) 

with relatively small wavenumbers, = 2n/kQ, and large width Qq, 
that is, Xq <ti Nr and Gq > Nr. Along an elastic chain (s = g = 0) 
with tn+m,n " 1 for m = ±1 and tn+m,n — 0 otherwise, the group 
velocity of the wave-packet^^'^^ is Vp = 2sin ko. 

Nonlinearity. The transmission coefficients Tlr (left- right) and Trl 
(right-left) have been computed as the ratio of the norm of the wave- 
packet transmitted, respectively, to the right and left chain, to the 
(conserved) total initial norm (see Methods for computational 
details). In Figure 2a Tlr and Trl have been plotted as functions 
of the magnetic flux, for different values of the V[\|/n] parameters a 
= (s,g) on the ring sites. Here, s is the same on both sides of the ring, 
5e = 0, whereas g is not, 5g = 0.3. The role of nonlinearity combined 
with asymmetry is apparent. Contrary to the symmetric and 
harmonic cases with go = 0 but SqT^O (dashed curves), the 
transmission coefficients across a nonlinear and upside-down 
asymmetric ring, Tlr(O) (solid, blue) and Trl(O) (solid, red), do 
depend on the direction of propagation. 

Qualitatively, this can be explained as an interference effect. 
Entering the ring from the left (right) at the junction site — Ni 
(+Ni), the wave-packet splits into two wave-packets (not necessarily 
with the same norm), which then propagate respectively along the 
upper and lower path. The magnetic flux changes their phase with 
opposite sign. As Qq > Nr, when the two wave -packets recombine at 
the exit junction sites Ni ( — Ni), they interfere. For a given ^ their 
phase difference at recombination depends on the parameters of the 
upper and lower ring sides. As long as the on-site potential is har- 
monic, go = 0, changing the direction of propagation only reverses 
the sign of the phase difference, so that Tlr(O) = Trl(O). For go^O, 
however, site resonances are detuned differently for ko > 0 and ko < 
0, which leads to different transmission coefficients^^. 

The O dependence of the transmission coefficients in the har- 
monic case, go = 0, is also remarkable. For So < 2 the curves 
Tlr,rl(^) are bimodal with symmetric maxima (close to 1) for So- 
dependent values of the flux. On increasing So the two transmission 
peaks tend to merge into a lower single peak centered at ^ = 0. Note 
that in any case the transmission coefficients are appreciably larger 
than 0.5 (incoherent propagation regime) only in the central flux 
range, -0.25 < ^/^o < 0.25. 

To further characterize the asymmetric transmission properties of 
our wave diode, in Figure 2b we also plotted the rectification factor 
Rf, 

^ TLR(ko)-TRL(-ko) 

as a function of Its overall magnitude strongly depends on the 
nonlinearity constant go. On keeping the asymmetry parameters 
fixed, 5e = 0 and 5g = 0.3, the amplitude of Rf goes through a 
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Figure 2 | (a) Dependence of the transmission coefficients, Tlr (solid blue) and Trl (solid red), on the magnetic flux O in units of Oq- The parameters Sq, 
and go are reported next to each curve and 5g = 0.3. The dashed curves represent Tlr = Trl as obtained for go = 0, 5g = 0 and different So- (b) Rectification 
factor, Rf vs. 0/(I)o, for different go (see legend). The simulation parameters are go = 1, Sq = 1.2, 5g = 0.3 and 5g = 0. In both panels the parameters 
of the wave-packet are Iq = 0.6, no = 250, = 25 and ko = 7i/2, see equation (3), and the ring is as in Figure 1. 



maximum for an optimal value of go, while Tlr and Trl are 
drastically lowered at larger go (not shown). Too strong a nonlinear 
coupling tends to suppress rectification. 

Finally, we notice that inverting the sign of go in Figure 2b does not 
appreciably affect the overall magnitude of the Rf oscillations. 
However, the curve for go = — 1 develops a strongly fluctuating O 
dependence, which hints at a chaotic behavior^^. This feature has 
been detected any time So and go had opposite sign, namely when 
the nonlinear on-site potential, V[\|/n], was bistable. Another 
example of rectification for negative nonlinearity constant is dis- 
played in Figure 3. 

Flux as control parameter. Figure 3 illustrates how the rectification 
power of the diode can be effectively controlled by tuning the 



magnetic field. Here we set O so as to maximize the difference 
between Tlr and Trl (i.e. Rf) and visualized the time evolution of 
the incoming wave-packet. When the packet proceeds from the left, it 
passes through the ring undergoing minor distortion; its reflected 
image is barely detectable. Vice versa, when the wave-packed impacts 
the ring from the right, it gets almost entirely reflected, only a small 
fraction of it being transmitted. During the transmission process the 
norm of the traveling wave tends to accumulate on the ring sites to 
different degrees. In the lower panels of Figure 3 we plotted the 
probability spectral density (p.s.d., defined in Methods) on site 
n = 0 for the two cases. The relevant p.s.d.'s peak respectively at 
zero frequency (wave not changing phase during the process), when 
the wave-packet gets transmitted, and at 27i;Vp/NR (wave shifting 
phase by n upon collision), when it gets reflected. A (nonlinear) 
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Figure 3 | Upper panels: Time evolution of the Gaussian wave-packet of Figure 2 injected respectively from the left (a) and right side (b) of the ring. 
The norm of the wave at each site of the device is represented according to the color code bar on the r.h.s. of the panels. Here, go = 1> ^/^o = ~0-2 and 
all remaining simulation parameters are as in Figure 2b. Lower panels: Corresponding probability spectral densities (p.s.d.) at site n = 0 in the same 
arbitrary units (see Methods). 
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Figure 4 | Diode performance on varying asymmetry and size of the AB 
ring: Rf vs. <l>/<I>o for (a) different 8g but 8^ = 0 and (b) different 6^ but 6g 
= 0; (c) Tlr vs. <I>/Oo for different ring circumferences Nr. In all panels Sq 
= 1.2, go = 1 and the injected wave-packet is as in Figure 2. Moreover, in 
(a) and (b) Nr = 8, and in (c) 5g = 0.3 and 5^ = 0. 

analysis of the fine structure of the peak in Figure 3d rests outside the 
purpose of the present report. 

Most remarkably, panels (a) and (b) of Figure 2 suggest that the 
symmetry relations Tlr(O) = Trl( — and, therefore, Rf{ — ^) = 
— Rf{C)) must hold for any choice of a. Indeed, under the transforma- 
tions a s^a, ko Skko, O S(d^, with Si = ± 1 and i = ot, k, ^, the 
dynamics of equation (1) with our choice of tn+m,n satisfies three 
distinct symmetry relations, which in compact form read (see 
Methods) 



Here, S(d = — 1 means that the sign of the AB phase (i.e., of the 
magnetic field) is reversed, 8k = — 1 indicates a change in the dir- 
ection of the incident wave-packet (i.e., an exchange of the suffixes 
LR and RL) and, finally, = — 1 corresponds to turning the on-site 
potential, V[\|/n], upside-down. Our numerical data confirm all iden- 
tities obtainable from equation (5). 

Asymmetry. In panels (a) and (b) of Figure 4 we analyze the role of 
the second key ingredient of our rectifier, namely its upside-down 
mirror asymmetry. In panel (a) 6g was increased at 5e = 0. The 
rectification factor in the domain —0.25 < ^/^o < 0-25 turned 
out to be rather sensitive to the asymmetry of the quartic term of 
V[v|/n]. The slope of Rf at O = 0 increases with 5g up to a maximum. 
At even larger 5g it changes sign and, as a consequence, a current 
inversion may occur as the curve Rf{0) oscillates around zero. 

Asymmetry can be controlled also by varying 5g at 5g = 0. As long 
as go ^ 0, this suffices in principle to operate our wave diode, as 
shown in Figure 4b. However, contrary to the situation discussed 
above, here increasing 5e induces frequent sign changes in the slope 
of Rf at 0 = 0, thus making an effective control of the diode operation 
more difficult. 

Finally, we briefly discuss the coherence issue. For the two split 
wave-packets in the AB ring to interfere at the exit junction site, 
where they recombine, is necessary that the relative spatial delay 
accumulated by following respectively the upper and lower path, 
be much smaller than their width. In Figure 4c we display Tlr as a 
function of the flux for different circumferences, Nr, of the AB ring. 
The O oscillations of the transmission coefficient get suppressed on 
increasing Nr and, eventually, for Nr > Gq Tlr fluctuates around its 
incoherent regime value, Tlr =1/2, which confirms the absence of 
interference between the two exiting wave packets. 

Discussion 

The wave diode we discussed in this report consists basically of an AB 
ring made of two halves with different nonlinear transport prop- 
erties. Its operation principle should not be mistaken for the standard 
rectification mechanism demonstrated in the recent literature^^, 
because it has especially been designed to retain longitudinal mirror 
symmetry. Here, it is the combined action of transverse mirror asym- 
metry with respect to the direction of transmission and the AB phase 
acquired by the charges moving around ring, that causes a net lon- 
gitudinal transmission current. The rectification power of the result- 
ing wave diode is controlled by the magnetic flux threading through 
the ring. The special topology of our AB diode could be implemented 
ad hoc, for instance, by intertwining two different molecular wires or 
unzipping a DNA chain so as to create a one dimensional molecular 
filament with interspersed loops. 

More interestingly, our wave-diode design suggests more direct 
applications to real electronic systems. Let us consider, for instance, a 
thin strip of a type II superconductor sandwiched between two strips 
of different conducting materials to form a three-layered long ribbon. 
On applying a constant magnetic field oriented parallel to the ribbon 
layers and perpendicular to the ribbon axis, one can generate one or 
more magnetic flux lines (fluxons), each of magnitude Oq, which 
under appropriate fabrication conditions, can be trapped at special 
pinning sites in the central layer. The intensity of the applied mag- 
netic field controls the fluxon density. Electrons moving in the outer 
ribbon layers are affected by the presence of such tightly sandwiched 
fluxons^^. Indeed, on passing each flux line, they acquire an AB phase, 
which depends in sign and magnitude on which conducting strip 
they propagate along. On assuming that the material of at least one of 
the two conducting strips has nonlinear transport properties, the 
situation envisaged in our stylized tunable Aharonov-Bohm diode 
is thus actually realized. 

Let us consider now an even simpler conducting ribbon by remov- 
ing the inner superconducting strip. At the interface between the two 
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remaining conducting strips, stresses induced, say, by mechanical or 
thermal deformation, can result in the nucleation of linear defects, 
dislocations, oriented perpendicular to the ribbon axis. Electrons 
moving past an edge dislocation accumulate an additional phase 
due to the distortion of the lattice structure that hosts the disloca- 
tion^^. In this situation the role of the magnetic flux lines would be 
played by the dislocations - both fluxons and dislocations can be 
regarded as different species of linear defects. In this setup, however, 
tunability would more difficult as dislocation density and topology 
are hardly controllable factors. A similar geometry could also work as 
phonon rectifiers provided tht the edge dislocations are replaced by 
screw dislocations and the incident packet is made of acoustic 
waves^^. 

Methods 

Simulations. In our simulations the number of ring sites was Nr = 4Ni, and the 
number chain sites on both sides of the ring Np = 500. The wave-packets travelling to 
the left (right) were injected at time t = 0 with width Gq = 25 and centered at no = 
Np/2 (-Np/2). At the endpoints of the right (left) chain, n = ±Np, we set tm+n,n = 0 
for m = ±1 (reflecting boundaries). For the purposes of the present study these 
chains can be considered infinitely long. Due to their large width the distortion of the 
transmitted (reflected) wave-packets was negligible (see, e.g.. Figure 3). 

The explicit expressions employed to compute the transmission coefficients were 



n = M, 
-N, 

n=-Np 



-N, 
n = N, 



(6) 



where the summations at the denominator have been computed at the injection time t 
= 0, and those at the numerator only after the time interval, (Np + Ni)/Vp, the 
transmitted wave-packet takes to reach the midpoint of chain on the opposite side of a 
ring with s = g = 0. In our simulations the residual norm trapped inside the ring was 
negligible for any choice of s± and g±. 

With regard to the denominators in equation (6), we notice that, while the nor- 
malization constant of the wave-packet at time t = 0, equation (3), has been set Iq = 
0.6 throughout our simulation work, the results reported in Figures 1-4 hold also for 
different normalization constants, Iq' = kIq, provided the nonlinearity constant go is 
rescaled as go go' = K^go. In particular, the wave-packet of equation (3) can be 
normalized to one incident particle by setting k = 0.15. 

The probability spectral densities (p.s.d) at site n = 0 were computed by means of 
the Fast Fourier Transform (FFT) utility of the software package OriginProS. The 
time series of the wave function \|/o(t) used in our computations were recorded 
between the time when the incident wave packet hit the injection ring site and the 
time when the center of the transmitted wave packet reached the midpoint of the exit 
chain. 

Model symmetries. The symmetry relations of equation (5) have been demonstrated 
numerically, but could well have been predicted on the basis of simple analytical 
arguments. 

1 . Rf( - O) = - Rf{0) follows from the identity Tlr(O) = Trl( - O), see Figure 2a, 
which is due to the right-left mirror symmetry of the device with respect to its 
center. Moreover, replacing n by -n in equation (1), however, inverts the sign of 
0). 

2. Rf( — ko) = — Rf(ko) follows from the very definition of Rf, equation (4), and the 
observation that by inverting the sign of ko one also interchanges the diode 
injection and exit site. 

3. Rf( — 8,— g) = — Rf{s,g) is a property of the cubic DNSE. Equation (1) describes 
the same transport dynamics even if one takes its complex conjugate after sub- 
stituting \|/n by \|/n*. Such a transformation would correspond to inverting the 
signs of Sn, gn> and taking the complex conjugate of tm+n,n- Taking the complex 
conjugate of the hopping constant can be restored by replacing ko with — ko. 
Numerical evidence of this symmetry relation is not reported in this work. 
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